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SUMMARY 

The traditional steady-state model of turbulent convection in a thermal boundary 
layer has been modified to include a nonsteady penetration component of heat transfer. 
The penetration mechanism is assumed to result from appreciable changes in the spec- 
ific volume of local agglomerates of fluid near the wall under heating conditions. In 
some respects the penetration mechanism is similar to boiling. The overall model has 
been applied to forced convection of several fluids near their critical thermodynamic 
state. Some success has been achieved in using the model to account for the differences 
between the experimental data and conventional turbulent flow heat-transfer correlations 
for variable property fluids. 


INTRODUCTION 

The peculiar forced- convection heat transfer by a fluid when it is near its thermody- 
namic critical point has been observed by numerous researchers (refs. 1 to 9). Needless 
to say, an explanation for these peculiar characteristics of heat transfer has been of con- 
siderable interest to the heat -transfer community. 

For the most part, explanations of this behavior have centered on the unusual phys- 
ical and transport properties of near- critical fluids. Efforts to correlate these data have 
involved manipulation of the fluid and transport properties as they are inserted into 
standard forced- convection correlations. In some cases, such as references 10 to 13, 
more elaborate analytical techniques were tried where variable property boundary layer 
analyses were applied. Such differential descriptions of the boundary layer, as contained 
in reference 10 or 14, were utilized. 

While limited success has been realized in correlating some near- critical fluid heat- 
transfer data by the selection of reference fluid and transport properties, general success 



by this method cannot be claimed. Consequently, other approaches need examination in 
rendering an explanation for the heat-transfer rates. 

One obviously different approach is to attribute the heat-transfer behavior to a change 
in the turbulent mechanism. In 1954, Goldmann suggested that a boiling-like mechanism 
of the near-critical fluid may explain the heat transport. Some visual studies described 
in references 15 to 17 have lent credance to such an explanation. The boiling-like explan- 
ation was extended further by Hsu in his written comments to reference 6. 

In 1966, Hall, Jackson, and Khan (ref. 18) suggested that the turbulent diffusion in 
the boundary layer may be amplified by the expansive coefficient of the fluid. They ar- 
gued that they could relate their carbon dioxide heat-transfer data to such a mechanism 
and not to the near- critical values of the transport properties. 

Herein, a somewhat different explanation which is based on a penetration model of 
the boundary layer is prof erred. Hanratty (ref. 19) first suggested that the turbulent 
boundary layer could be represented by a penetration model. In the conventional turbu- 
lent heat -transfer approaches, the fluctuating components of velocity and temperature 
are represented by some statistical time-averaged value. Through averaging, one can 
eliminate time as a variable. In the penetration model approach, however, the time de- 
pendence is emphasized. The heat transfer is treated as a transient conduction process, 
which is repetitive as fluid species contact the wall. As is well known, the penetration 
model had its origin in gas-liquid diffusion processes and has been widely used in compu- 
tations of processes in fluidized and liquid beds. Its more recent application to the tur- 
bulent boundary layer has been justified by evidence of turbulent -like behavior in the so- 
called laminar sublayer. Many references, such as 20 to 24, contain evidence of this 
sublayer turbulence which has been labeled by some as a "sublayer instability. " 

The heat-transfer model incorporating a penetration component is compared to the 
liquid hydrogen experimental data found in reference 25. These near-critical data have 
not been adequately correlated by standard correlation schemes. However, the relative 
success in correlating these data is used as a supporting argument for the introduction 
of a new mechanistic effect. The correlation equation developed for the hydrogen data 
is applied to a limited amount of heat-transfer data for near-critical carbon dioxide and 
water. The main purpose of this report is to present a proposal which takes into account 
the peculiar turbulent behavior of near- critical fluids. 


SYMBOLS 


A area 

a frequency enhancement factor 

C constant 

2 



c specific heat 

D diameter of channel 

d diameter of jet (appendix B) 

G mass flow per unit area 

g gravity 

H enthalpy 

h heat-transfer coefficient 

J mechanical equivalent of heat 

k thermal conductivity 

L length 

Nu Nusselt number 

n exponent 

Q heat rate 

q heat flux 

Re Reynolds number 

St Stanton number 

T temperature 

U velocity component, axial 

V velocity component, normal 

v specific volume 

x axial distance and quality fraction 

y normal distance 

p density 

co frequency 

p viscosity 

cp thermal boundary layer thickness 

t time 

Subscripts: 

av average 

b bulk 



c contact 

enh enhanced 

exp experimental 

f film 

fg film reference based on gas 

fm film reference based on two-phase mixture 

i wall temperature, inside pipe (eq. (6)) 

o reference, adiabatic 

p penetration model 

t total 


PROPOSED HEAT-TRANSFER MODEL 

In view of the enhanced heat -transfer rate observed with near-critical fluids (as com- 
pared to a film correlation), it is proposed that this enhancement be represented by an 
additive term to the conventional film property forced- convection coefficient. In general, 
the average heat-transfer coefficient can be written as 


h av ~ C l h f + C 2 h enh 


( 1 ) 


One can introduce an enhanced coefficient and the appropriate weighting coefficients Cj 
and Cg to enable correlation of the experimental data. 

In this report it is proposed that the enhancement mechanism is "penetration.” It 
is analogous to the penetration mechanism of packed beds in chemical engineering tech- 
nology. The extreme density gradient through the boundary layer produces instabilities 
which disrupt the viscous, near-wall portion of the turbulent boundary layer. The pene- 
tration of this near-wall sublayer by fluid packets from the outer region of the boundary 
layer is the reason for the name of the mechanism. The weighting coefficients reflect 
area fractions of the surface where the penetration mechanism does and does not domin- 
ate. 


A gross overall view of the penetration model and the involvement of an area- 
weighting factor is shown in figure 1. Equation (1) can be written particularly for the 
penetration enhancement model as follows: 


av 



(2) 


4 


Dissipation 




Figure L - Concept of simultaneous hydrodynamic and penetration mechanisms. 


where 

Ap/Af penetration area fraction 

h penetration coefficient 

r 

conventional forced convection coefficient, evaluated at film properties 

In this approach, the area fraction parameter and the penetration heat -transfer coef- 
ficient are only correlation parameters. The magnitude of these coefficients was estab- 
lished by the use of a limited amount of near-critical hydrogen data from references 25 
and 26. Data points from 11 runs were used in determining the constants. The correla- 
tion method was eventually applied to a total of 50 hydrogen data runs; data from several 
axial stations were incorporated for each run. The run numbers for the 11 runs used in 
determining the constants are marked with an asterisk in tables I and II. The details of 
how the coefficients were estimated appears in appendix A. 
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In the prediction of the penetration heat-transfer coefficient, the essential elements 
are the bulk properties of the fluid and the frequency with which the packets disrupt the 
sublayer. The form of the equation is exactly similar to the one used to describe pene- 
tration heat-transfer coefficients in packed beds (ref. 30): 



where u> is the frequency and the subscript b denotes bulk conditions. 

As described in appendix A, the frequency is presumed dependent on density change 
in the boundary layer and on the average momentum of the axial flow in the boundary 
layer. The interaction of density change and flow momentum is considered in an analo- 
gous fashion to forced convection boiling. The behavior of the frequency term was com- 
puted from 


w = „ R p b - <v) V 2J g cK - H b) (4) 

0 L p f U » p b_ 

The term in brackets, involving the density change, axial momentum of the boundary 
layer, etc. , may be interpreted to be a dimensionless frequency enhancement factor. 

The term w Q in equation (4) is strictly an empirical frequency but it does turn out that 
its magnitude (20 cps) is within the range of frequency measurements made in the sub- 
layer under adiabatic conditions. The enhancement factor is always much greater than 
one, so the frequency computed from equation (4) always exceeds the reference value 
<x) Q . Equation (4) is only applicable to the heating case where p^ > p w< 

Hydrodynamic arguments about mass transfer have served to rationalize the corre- 
lation form for the penetration area fraction used in estimating the average heat-transfer 
coefficient. The penetration area fraction weighting term is primarily a function of the 
heat-transfer rate. The rate can be represented by the bulk Stanton number. The corre- 
lating equation for the penetration area fraction is 



This correlation scheme (described in greater detail in appendix A) was applied to the 
near-critical heat -transfer data for liquid hydrogen found in reference 25. As mentioned 
previously, 11 runs from reference 25 were utilized in determining the coefficients em- 
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pirically. Then, the correlation method was extended to most of the near-critical hydro- 
gen heat-transfer data found in reference 25. It is important to observe that an upper 
limit to the magnitude of Ap/ was assigned. The value was approximately 0. 16 and 
it is akin to the maximum area fraction observed on boiling surfaces near burnout con- 
ditions . 


COMPARISON OF MODELWITH EXPERIMENTAL DATA 

Hydrogen 

The mechanistic model of turbulent heat transfer as set forth in the previous section 
and appendix A was compared with experimental near-critical data for hydrogen, carbon 
dioxide, and water. The hydrogen data from reference 25 were the chief source of data. 
The empirical constants in the model were obtained from comparisons with selected 
near- critical hydrogen heat-transfer data from reference 25. 

It was stated in the INTRODUCTION that the near-critical forced-convection data for 
hydrogen could not be correlated by a standard turbulent correlation. As an example of 
the disparity between a standard correlation and experimental data, the heat-transfer 
coefficient distributions along the axis of a tube are shown in figure 2. The range of 
L/D spans from 13 to 34. The circled points are experimental data, and the squares are 
predicted values from a turbulent pipe flow heat-transfer correlation. Note the wide 
discrepancy between prediction and experiment for this run. 



Figure 2. - Comparison of film penetration model and experimental 
heat-transfer coefficients. 
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Also shown in the figure is the heat-transfer coefficient prediction employing the 
penetration model concept (eq. (2)). For this particular run, the prediction agrees well 
with experiment. Incidently, this run (3-1094) was not one of the original 11 runs used 
in determining the magnitude of the empirical constants. It does happen that the bulk 
fluid conditions are very close to critical for both temperature and pressure. 

As a check on the versatility of the model to fluid state conditions far removed from 
the critical, high-pressure hydrogen data involving appreciable wall- to bulk-temperature 
ratios were included in the comparisons. Selected runs from references 26 and 27 were 
used for that purpose. Tables I and II are a listing of the run numbers and test section 
stations from references 25 to 27 and a tabulated comparison of the heat-transfer coeffi- 
cients with the experimental values for hydrogen. 

A quick observation of the series of graphs in figure 3 shows the relative merits of 


Pressure, 

psia (N/cm^) 

O 250(or less) (172) 

A 275 (189) 

O 300 (206) 

□ 400 (275) 


"8 

E 
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(a) Large diameter tubes, d, 0. 438 and 0, 507 inch ( 1. 11 and 1. 29 cm). 



O 


(b) Medium diameter tubes, d, 0.335 inch (0.85 cm). 
3[- 


2\— O 



Ratio of experimental to film correlation 
heat-transfer coefficients, h exp /hf 

(c) Small diameter tubes, d, 0. 220 inch (0. 56 cm). 


Figure 3. - Comparison of penetration model and film correlation 
predictions to experiment for cryogenic hydrogen. 
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TABLE I. - COMPARISON OF MODEL TO HYDROGEN EXPERIMENTAL INFORMATION 





































































TABLE I. - Continued. COMPARISON OF MODEL TO HYDROGEN EXPERIMENTAL INFORMATION 




(b) Test section diameter, 0. 335 inch (0. 85 cm) (ref. 25) 


Run 

Axial dis- 
tance to 
diameter 
ratio , 
L/D 

Bulk 

pressure, 

p b ’ 

psi 

Wall 

temperature, 

T , 
w’ 

°R 

Bulk 

temperature, 

T b- 

°R 

Heat flux, 

q, 

Btu/(in. 2 )(sec) 

Bulk axial 
velocity 
component, 

u b 

Area frac- 
tion of 
penetration, 

V A t 

Heat -transfer 
coefficient ratio 

Penetration 

contribution 

fraction, 

Vav/ \ t/ 

h exp /h f 

h exp /h av 

1-880 

19. 5 

250 

697 

64.6 

2.39 

105.6 

0. 100 

3.28 

1.17 

0.67 


28.5 



960 

70.0 



175.6 

.103 

2.35 

1.06 

.59 


34. 5 



906 

77.9 



240.2 

. 127 

1.86 

1.41 

.33 


46. 5 




835 

105 



376.9 

.158 

1.38 

1.40 

.13 

14-702 

19.5 

250 

465 

62.5 

1.43 

66.5 

0.115 

3.45 

1.32 

0.70 


28. 5 



672 

64.7 



99.8 

. 113 

2.42 

.87 

.67 


34.5 



757 

67.7 



133.4 

.118 

1.98 

.90 

.59 


46.5 



709 

80. 1 



216 

.158 

1.39 

1.18 

.26 

16-700 

19. 5 

275 

432 

59.1 

1. 50 

84 

0.078 

3.00 

1. 55 

0.52 


28. 5 



572 

62.4 



100.2 

.084 

2.47 

1.09 

.59 


34. 5 



670 

63.8 



117.4 

.085 

2.18 

.95 

. 59 


46. 5 



597 

66. 2 



167 

. 122 

1.68 

1.02 

.45 

*30-886 

19. 5 

300 

552 

67.4 

1.87 

79.2 

0. 120 

3.02 

1.02 

0.70 


28. 5 



773 

73.2 



121.6 

.121 

2.21 

1.10 

.52 


34. 5 



870 

80.9 



161.6 

. 129 

1.83 

1.20 

.42 


46. 5 



854 

103. 5 



252 

.156 

1.36 

1.26 

.20 

*55-258 

19. 5 

400 

383 

59.4 

0.67 

44.5 

0.072 

1.91 

1.00 

0.51 


28. 5 



478 

63. 7 



49.8 

.080 

1. 65 

.71 

.60 


34. 5 



458 

66.1 



54.4 

.098 

1.57 

.67 

.61 


46. 5 



396 

69.8 



65. 8 

. 146 

1.44 

.71 

.56 

57-881 

19. 5 

700 

952 

77.9 

2.23 

82.0 

0.063 

1.99 

0.94 

0.55 


28. 5 



887 

87.9 



105,9 

.093 

1.64 

1.05 

.41 


34. 5 



819 

94.9 



124,8 

.117 

1.45 

1.06 

.34 


46. 5 



731 

110.4 



165.3 

.156 

1.20 

1.09 

.22 

4-888 

19. 5 

250 

898 

68.7 

2.23 

120.9 

0. 104 

1 

2.85 

1.12 

0.64 


28. 5 



1099 

84.7 



202.2 

. 120 

2.01 

1.25 

.44 


34. 5 



1053 

102.4 



271.5 

. 133 

1.64 

j 1.33 

.29 


46. 5 



1003 

139.8 



402.7 

.158 

1.29 

1.29 

.15 

3-887 

19. 5 

250 

628 

63.2 

2. 

51 

111.9 

0.088 

3. 32 

1.39 

0.62 


28. 5 



889 

65.7 



166.2 

.088 

2. 57 

1.13 

.55 


34. 5 



862 

68.7 



222.5 

.106 

2.06 

1.29 

.43 


46. 5 

\ 


737 

79.9 1 

i 


354.2 

. 158 

1.46 

1. 39 

.16 

43-688 

19. 5 

500 

692 

68.4 

1. 

22 

54.7 

0.063 

2.01 

0.84 

0.61 


28. 5 



771 

74.6 



70.0 

.078 

1.72 

.71 

.62 


34. 3 



662 

78.6 j 



87.6 

. 10 

1.49 

.77 

.53 


46. 3 



550 

87. 5 



111 

. 14 

1.24 

.94 

.35 

47-653 

19. 5 

400 

694 

60.2 

1. 

3 

86.1 

0.04 

1.79 

1.11 

0.40 


28.4 



642 

64.7 



97.8 

.06 | 

1.64 j 

. 95 

.45 


34.3 



553 

67.0 



107.2 

.075 

1.56 

.95 

.44 


46.3 



464 

70. 5 



130.3 

.10 

1.40 

1.00 

.36 

32-312 

19.4 

325 

455 

63.9 

1. 

82 

75 

0.098 

3.37 

1.44 

0.61 


28.35 



685 

67.7 



104 

.096 

2.48 

.97 

.64 


34.3 



759 

75.2 



171 

.118 

1.75 

1.15 

.42 


40 



666 

88. 5 



252 

.135 

1.32 

1.25 

.18 

34-273 

19.4 

400 

595 

63.1 

1. 

82 

94.7 

0.058 

2.38 

1.45 

0.42 


28.3 



585 

69.3 



124 

.087 

1.97 

1.20 

.44 


34.3 



567 

70.5 



133.6 

.097 

1.87 

1.21 

.42 


40 



473 | 

75.4 



203.8 

. 131 

1.58 

1.32 

.27 


*Run numbers used in determining constants. 
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TABLE I. - Concluded. COMPARISON OF MODEL TO HYDROGEN EXPERIMENTAL INFORMATION 


(c) Test section diameter, 0.22 inch (0. 56 cm) (ref. 25) 


Run 

Axial dis- 

Bulk 

Wall 

Bulk 

Heat flux, 

Bulk axial 

Area frac- 

Heat -transfer 

Penetration 


tance to 

pressure, 

temperature, 

temperature, 

q , 

Btu/(in. 2 )(sec) 

velocity 

tion of 

coefficient ratio 

contribution 


diameter 




component, 

penetration, 



fraction, 



1 w ’ 

i b’ 

Wf 

h exp/ h av 


ratio, 

psi 

°R 

°R 


u b ' 

VA 

(\V a p\ 


L/D 





ft/sec 



Uw 

3-1022 

29 

300 

748 

67.1 

2. 87 

283 

0.076 

1.78 

1.38 

0.28 


36 

| 

690 

69.3 

| 

355.7 

.097 

1. 55 

1.39 

.19 


45 

l 

601 

73.8 

* 

473.1 

.135 

1.36 

1.39 

.09 

9-1043 

29 

225 

201 

56.9 

0.76 

141.1 

0.083 

2.42 

2.11 

0.19 


36 

1 

265 

57,9 

1 

146.7 

.077 

1.88 

1.54 

.23 


45 

4 

315 

59.1 

1 

155.3 

.080 

1.55 

1.23 

.25 

11-1026 

29 

300 

775 

78.2 

2.87 

231 

0.09 

1.63 

1.36 

0.26 


36 

I 

732 

84 

I 

286 

.117 

1.42 

1.34 

.18 


45 

♦ 

660 

94 


367 

.156 

1.24 

1.32 

.10 

*22-1029 

29 

700 

851 

83.7 

2.88 

147 

0.069 

1.48 

1.20 

0.24 


36 

I 

712 

93.4 

| 

187 

.111 

1.25 

1.18 

.16 


45 


667 

100. 5 

♦ 

217 

.139 

1.21 

1.19 

.14 

5-856 

29 

400 

783 

70.7 

3.3 

196 

0.068 

1.96 

1.45 

0.31 


34 

1 

688 

76.2 

1 

269 

. 142 

1.60 

1.41 

.20 


45 

1 

635 

83.3 

1 

363 

.156 

1.34 

1.36 

.10 

30-1038 

23 

700 

178 

61.2 

1.39 

150.2 

0.109 

1.56 

1.67 

0.03 


33 

1 

197 

65.1 

1 

158.9 

.137 

1.49 

1.62 

.10 


40 

4 

191 

67.4 


165.4 

. 158 

1.52 

1.69 

.06 

20-1036 

23 

» 

187 

59.3 

1.35 

184 

0.088 

2.12 

2.18 

0.06 


33 

mm 

244 

61.8 

1 

196 

.095 

KIM 

1.73 

.10 


40 

mm 

226 

63.3 

i 

205 

.12 


1.84 

.10 

4-1007 

29 


490 

62.2 

WBM 

, ■ 

0.064 

1.62 

mm 

0.28 


36 


481 

63.1 

warn 

1 

.077 

1.53 

m 

.28 

2-1005 

29 


774 

68.9 

2.71 


0.082 

1.83 

1.37 

0.31 


36 


708 

70.6 

1 

286 

.095 

K£i 

1.42 

.19 


45 

■■ 

630 

80. 5 

1 

439 

.13 

Hi 

1.39 

.09 

14-1025 

29 

400 

722 

72 

2.88 

219 

0.073 

1.67 

1.32 

0.27 


36 

1 

679 

75. 1 

1 

262 

.091 

1.50 

1.30 

.21 


40 


594 

79.8 

4 

330 

.11 

1.35 

1.33 

.12 


*Run numbers used in determining constants. 


the model in predicting convective heat transfer for near-critical hydrogen as compared 
to the conventional forced convection correlation. The three graphs shown segregate the 
results according to tube diameter - large, medium, and small. This is the same as 
segregating the data according to ranges of free stream velocity. The magnitudes of 
free stream velocities were inversely proportional to the test section diameters. The 
ordinates are ratios of the experimental heat -transfer coefficients to the predicted coef- 
ficients from the model. The abscissas of the figures are the ratios of the experimental 
heat-transfer coefficients to the coefficients predicted by the conventional Dittus-Boelter 
type correlation. The most dramatic results are associated with the large tube diameter 
data. Extremely large ratios of experimental heat -transfer coefficient to film correla- 
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TABLE n. - HIGH-PRESSURE HEAT-TRANSFER DATA 


(a) Test section diameter, d = 0. 21 inch (0. 53 cm) (ref. 27) 


Run 

Axial dis- 
tance to 
diameter 
ratio , 
L/D 

Bulk 

pressure, 

P b’ 

psi 

Wall 

temperature, 

T w- 

°R 

Bulk 

temperature, 

T b’ 

°R 

Heat flux, 

q, 

Btu/(in.^)(sec) 

Bulk axial 
velocity 
component, 

u b> 

ft/sec 

Area frac- 
tion of 
penetration, 

y\ 

Heat-transfer 
coefficient ratio 

Penetration 

contribution 

factor, 

/h \/A \ 

1 Ml P 1 

\ h av/W 

h exp/ h f 

h exp //h av 

5 

9.5 

1000 

1326 

82.7 

12.4 

291 

0.03 

2.30 

2.10 

0.12 


19 



1296 

100.1 

11.9 

394 

.06 

1.82 

1.72 

.11 


28.6 



1030 

118 

10.98 

521 

.09 

1. 55 

1.60 

.07 


38 



855 

137 

10.7 

654 

.14 

1 . 45 

1.62 

.04 

6 

9.5 

1000 

864 

74.6 

9.7 

279 

0.04 

2.19 

2.15 

0.05 


19 



1253 

87.8 

10.4 

346 

.05 

1.76 

1.62 

.12 


28.6 



980 

100.6 

9.4 

438 

.08 

1.53 

1.52 

.09 


38 



773 

113.5 

9.0 

544 

.12 

1.42 

1.54 

.04 

28 

9.5 

1250 

1093 

81.7 

10.2 

208 

0.04 

2.18 

2.01 

0.11 


19 



1370 

100.6 

10.9 

265 

.06 

1.80 

1.60 

.16 


28.6 



806 

119.2 

10.2 

337 

.13 

1.85 

1.97 

.07 


38 



847 

137 

9.3 

412 

.06 

1.47 

1.63 

.06 

35 

9.5 

550 

1378 

70.3 

7.4 

224 

0.03 

2.46 

1.83 

0.27 


19 



954 

79 

6.6 

299 

.06 

2.11 

1.78 

.20 


28.6 



755 

87.7 

6.2 

405 

.09 

1.76 

1.77 

.09 


38 



663 

97.7 

5.7 

527 

.13 

1.44 

1. 55 

.05 


(b) Test section diameter, d = 0.22 inch (0. 56 cm) (ref. 26) 


206-A-6 

21 

26 

30.8 

1450 

I 

930 

904 

961 

92 

99 

106 

5.6 

i 

189 

205 

221 

0.06 

.08 

.08 

1.31 

1.26 

1.18 

1.22 

1.19 

1.11 

0.13 

l 

103-A-4 

14.3 

750 

957 

143 

1.9 

121 

0.07 

1.23 

1.04 

0.21 


17.5 

750 

1098 

151 

1.9 

130 

.07 

1.11 

.92 

.23 

230-A-6 

21 

1500 

708 

126 

5.6 

221 

0.08 

1.21 

1.27 

0.04 


28 

1500 

659 

133 

5.6 

281 

.10 

1.18 

1.26 

.04 

♦388-A-7 

17.1 

1350 

442 

87 

0.9 

34 

0.09 

1.28 

0.94 

0.33 


19.4 

1350 

471 

89 

.9 

35 

.09 

1.21 

.86 

.35 


(c) Test section diameter, d = 0.335 inch (0.85 cm) (ref. 26) 


171-A-4 

20.6 

900 

323 

67 

2.8 

209 

0.07 

1.47 

1.53 

0.03 


30.2 

| 

388 

72 

I 

225 

.08 

1.24 

1.28 

.05 


33.4 

1 

351 

74 

1 

231 

.10 

1.31 

1.37 

.05 


♦Run numbers used in determining constants. 


tion coefficient were observed. The relatively small magnitude of velocity in the large 
tubes probably enabled the penetration mechanism to go on relatively unabetted by free 
stream momentum effects. If one judges from figure 3(a) (large diameter tubes), it does 
appear that the penetration model concept is quite successful in reducing the differences 
between experimental and predicted values of the heat-transfer coefficient. Similar, but 
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lesser, improvements in the prediction of the heat-transfer coefficients were observed 
with smaller diameter tubes (figs. 3(b) and (c)). 

An appreciation of the relative magnitude of the penetration mechanism contribution 
to the overall heat -transfer mechanism can be discerned from figure 4, which compares 
the ratio of the heat flux contribution ascribed to penetration to the average overall heat 
flux with the ratio of the experimental heat-transfer coefficient to the film coefficient. 
Over a considerable range of conditions, the relative contribution of the penetration 
mechanism does increase as the ratio of the experimental heat-transfer coefficient to the 
film coefficient increases. The relative contribution of the penetration mechanism ap- 
pears to peak at approximately 80 percent of the predicted average heat-transfer rate. 
Note that the near-critical pressure data make up the high end of the plot. 
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Figure 4. - Relative contribution of penetration mechanism to overall 
heat transfer. 


The reader may wonder if this figure just constitutes a sort of back- calculation dis- 
play of some of the empirical aspects of the model. It is true that the magnitudes of the 
coefficient in equation (A7) (expression for hp) and Cg in equation (A9) (expression for 
Ap/A t ) were determined from comparisons with a few sample runs. After fixing these 
constants, equations (A7) and (A9) were applied to a large array of hydrogen data over a 
considerable range of conditions. Computed parameters, such as the penetration coeffi- 
cient hp, varied over one order of magnitude. Thus, the author claims that the method 
applies to a broad range of conditions where density effects, primarily, alter the conven- 
tional representation of turbulent forced convection. 
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The correlation equation was further compared to some high pressure, high heat 
flux data for liquid hydrogen presented in references 26 and 27. The pressure of the 
liquid hydrogen was well above critical, but the temperature was around critical. For 
these conditions, gradients in density and transport properties are not as severe as near 
the critical point. 

Table II is a tabulation of the predicted heat-transfer coefficients compared to the 
experimental values of references 26 and 27. Again the conventional film correlation 
and the penetration model correlations are presented. Cne can conclude for these data 
that the film correlation and the correlation involving the penetration model predict about 
the same values. The reason for this is that the penetration component for these super- 
critical pressures is small. Thus, the "average” coefficient is principally the film coef- 
ficient (see eq. (2)). 


Carbon Dioxide and Water 

A limited comparison of the model predictions with experimentally determined heat- 
transfer data was made for carbon dioxide and water. Data from references 8 and 9 
were utilized. The information on the run numbers and comparisons of the predictions to 
experimental values are presented in table III. The empirical constants derived from the 
hydrogen calculations were employed in the computational procedures for these fluids. 

Regarding the carbon dioxide data, the experimental heat fluxes are very low and the 
ratios of wall to bulk temperatures are also very low. The experimental data points were 
obtained at one station in the test section. The runs selected included data with bulk tem- 
peratures above and below the critical value. 

Actually, the operating pressure was slightly above the critical pressure, so the 
bulk temperature range represents conditions below the transposed critical temperature 
(locii of maximum C ) and at the transposed critical temperature. For runs 23 and 25 
the bulk temperature was less than the critical value and there appears to be no enhance- 
ment of the heat transfer. When experimental conditions were altered such that the bulk 
temperature was essentially at the transposed critical value, the heat-transfer coefficient 
reached a maximum value. However, the predicted film correlation coefficient was far 
below the experimental value. The introduction of a penetration component improved the 
prediction considerably as is evident from an examination of runs 28 and 30. 

About all that can be concluded from this comparison with experimental carbon 
dioxide heat-transfer data is that several experimental data points were predicted fairly 
well by the analysis. For these cases the predictions were generally superior to esti- 
mates obtained from pipe flow correlations. Also, the empirical constants of the model, 
which were derived from hydrogen data, appear to be applicable to another fluid. 
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TABLE HI. - COMPARISON OF MODEL TO EXPERIMENTAL INFORMATION FOR CARBON DIOXIDE AND WATER 


(a) Fluid, carbon dioxide; test section diameter, d = 0.9 inch (2.3 cm) (ref. 8) 


Run 

Axial dis- 
tance to 
diameter 
ratio, 
L/D 

Bulk 

pressure, 

P b’ 

psi 

Wall 

temperature, 

T w> 

°F 

Bulk 

temperature, 

T b' 

°F 

Keat flux, 

q, 

Btu/(in. 2 )(sec) 

Bulk axial 
velocity 
component , 

u b- 

ft/sec 

Area frac- 
tion of 
penetration, 

V A t 

Heat -transfer 
coefficient ratio 

Penetration 

contribution 

fractor, 

h exp/ h f 

^exp/ ^av 

23 

30.8 

1075 

92.9 

83.3 

0.072 

13.7 

0.156 

0.66 

0.71 

0.45 

25 





94.0 

85 



13.7 



1.61 

1.60 

.11 

28 





98.7 

87 



13.7 



3.70 

2.55 

.39 

30 





133 

88.8 



16.4 



2.14 

.91 

.64 


(b) Fluid, water; test section diameter, d = 0.371 inch (0.94 cm) (ref. 9) 


23200 

40 

— 

4000 

887 

731 

1.05 

16.9 

0.156 

2.37 

0.74 

0.74 

23202 

73 

4000 

922 

‘738 

1.05 

18.9 

.156 

2.23 

.73 

.72 

23206 

137 

4000 

904 

749 

1.06 

25 

.156 

2.04 

1.19 

.51 

22108 

40 

3300 

832 

691 

.57 

7.6 

.14 

3.31 

.42 

.89 

22112 

120 

3300 

820 

704 

.57 

10 

.14 

3.24 

.58 

.85 

22114 

137 

3300 

810 

707 

.58 

11.4 

.156 

3.09 

.61 

.83 

22445 

57 

3300 

923 

702 

0.37 

4.0 

0.156 

3.56 

0.10 

0.97 

22453 

186 



1131 

725 

.39 

10.6 

.156 

1.33 

.26 

.83 

22461 

120 



731 

700 

1.31 

3.6 

.156 

1.72 

.38 

.80 

22431 

25 



932 

687 

.76 

7.4 

.10 

3.66 

.38 

.90 

22134 

73 



1111 

703 

.77 

9.6 

.145 

2.67 

.16 

.95 

22874 

8 



746 

719 

.15 

41 

.03 

2.29 

2.1 

.08 

22875 

24 

3300 

747 

720 

0.15 

41 

0.07 

2.22 

1.9 

0.19 

22876 

40 

3300 

752 

721 

.15 

41 

.10 

1.99 

1.65 

.25 

22877 

57 

3300 

755 

722 

.15 

41 

.11 

1.89 

1.55 

.26 

23222 

8 

4000 

787 

747 

.24 

27 

.032 

2.39 

1.99 

.19 

23223 

24 

4000 

791 

748 

.24 

28 

.076 

2.25 

1.51 

.38 

23224 

40 

4000 

796 

749 

.24 

28 

.11 

2.07 

1.21 

.48 

23225 

57 

4000 

800 

751 

0.24 

29 

0.11 

1.95 

1.14 

0.48 

22898 

8 

3300 

806 

720 

.34 

35 

.04 

4. 57 

2.70 

.43 

22899 

25 



819 

723 



40 

.08 

4.13 

2.0 

.55 

22900 

40 



834 

726 



45 

.12 

3.70 

1.79 

.57 

22901 

57 



846 

730 



50 

.13 

3.43 

1.96 

.50 


The near- critical water heat-transfer data from reference 9 were utilized in a com- 
parison. The authors of reference 9 correlated their data by a regression technique and 
suggested the following correlation equation: 



0.00459 




Mi 

k i 


0.613 



( 6 ) 


where H is enthalpy, v is specific volume, and the subscripts i and b refer to inside 
wall conditions and bulk conditions, respectively. They claimed that this form of convec- 
tive correlation enabled better comprehension of the widely varying physical and trans- 
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port properties than the traditional pipe flow correlation wherein properties are evalu- 
ated at some reference condition. 

It is interesting to note that the near- critical water heat-transfer data exhibit the 
same trends as the hydrogen and carbon dioxide data: that is, the experimental heat- 
transfer coefficient is always greater than the predicted value based on film properties. 

Several runs from reference 9 are tabulated in table 111(b) and the data are compared 
to the penetration contribution model. For the low velocity runs (u < 10 fps), the pene- 
tration model overpredicted severely. For the higher flow velocities, the penetration 
model underpredicted, but generally the degree of underprediction was less than what 
the film correlation provided. 

It cannot be claimed that the penetration model presented herein is a satisfactory 
means of predicting the near-critical water heat-transfer data of reference 9. Neverthe- 
less, the trends of enhancement do seem to be in agreement with the penetration explan- 
ation. 


APPLICATION OF PENETRATION MODEL 
TO TWO-PHASE HEAT TRANSFER 


One two-phase flow regime which appears applicable to the penetration model is 
mist flow. Tiny liquid droplets migrate transversely from the bulk to the wall where 
they appear to participate in conduction and evaporation mechanisms. These migrating 
droplets can be considered to be penetration packets. 

The heuristic model which was developed for the near- critical fluid was applied to 
experimental mist flow data. The same equations and constants were utilized except for 
the area fraction parameter. The area fraction was calculated from estimates of the 
liquid void fraction computed in reference 28. Also, a film density was utilized which 
was weighted with the mass fraction of liquid and vapor at film conditions. This refer- 
ence density labeled p^ was defined to be 


p fm ~ 


1 

x 1 - x 


p fg p b 


(7) 


where x is the estimated average quality, p ^ is gas density evaluated at film temper- 
ature, and p b is liquid density. 

With the exception of these modifications, the approach employed for the near- 
critical fluids was applied to the two-phase problem. The correlation coefficients were 
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identical when equations (A3) and (A4) were employed. Figure 5 compares the estimate 
involving the penetration model with the film correlation approach. It is apparent that 
the penetration model contribution improves the comparison with experimental results. 



Ratio of experimental to film-correlation 
heat-transfer coefficients, h eX p/hf m 

Figure 5. - Comparison of penetration model for two- 
phase hydrogen to single-phase correlation. 
Diameter, d, 0.375 inch (0.95 cm). 


CONCLUSIONS 

Not withstanding the empirical nature of the description of the penetration mechan- 
ism, it does appear that the inclusion of such a component in a convective heat-transfer 
mechanism helps to explain the dramatic changes in heat transfer observed with near- 
critical fluids. It is further concluded that the peculiar heat transfer of this regime 
appears to be associated with the nature of the turbulence structure of the near-wall 
region. Fluid properties are important in the penetration mechanism. In particular, 
density and the rate of change of density of expansion mechanistically determine the fre- 
quency of the penetration cycle. Thus, a sizeable portion of the influence on overall heat 
transfer by fluid properties is carried into the analysis as a mechanistic effect. 

Recognition of the time-unsteady contribution to the heat transfer does emphasize 
that some comprehension of the instantaneous heat-transfer processes be achieved before 
any statistical averaging be done. The importance of assessing various mechanistic 
effects before averaging was also observed in studies of nucleate boiling (ref. 29). In the 
development of the standard turbulent flow heat-transfer correlations, time-averaged 
values of velocity and temperature are introduced initially. The mechanism is treated 
as a time-steady phenomenon which utilizes average values to comprehend the turbulent 
behavior. 

The model presented herein suggests that the heat-transfer mechanisms of the near- 
critical fluid and two-phase fluid (boiling) are similar in some respects. Such a similar- 
ity was suggested several years ago, and it has had experimental verification through 
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visualization techniques. It does seem that with further refinements a unified model of 
turbulent convection is possible for all the fluid states. This would obviate the current 
practice of abruptly changing heat-transfer models when the boundaries of fluid states 
are crossed. 

Finally, the model portrayed herein is a first approximation to a mechanism that 
has been qualitatively discerned experimentally. It is presented as one heuristic model 
to explain an enhanced component of heat transfer. Detailed information on the sublayer 
structure is needed to produce better quantitative data for verification and refinement of 
the model. The comparison of experimental data with a predicted heat-transfer coeffi- 
cient, which included a penetration contribution, was presented to support the arguments 
for the model proposed and not to propose another heat-transfer correlation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 29, 1969, 

129-01. 
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APPENDIX A 


HEURISTIC MODEL OF TURBULENT TRANSPORT WITH PENETRATION COMPONENT 

Penetration Heat-Transfer Coefficient 

The key concept in formulating the model of the penetration component of the bound- 
ary layer is borrowed from the literature dealing with the heat transfer in fluidized beds. 
The picture of a discrete fluid packet migrating toward the wall, as shown in figure 1, is 
similar to the model used in reference 30, for instance, to estimate the heat transfer 
between a fluidized bed and the vessel contact surfaces. As pointed out in the INTRO- 
DUCTION, Hanratty (ref. 19) must be credited with suggesting that such a penetration 
concept, which was originally applied to mass diffusion, could be applied to the turbulent 
boundary layer. 

For the estimation of hp, it is assumed that the packet penetrating to the wall 
carries its bulk properties and then begins its transient conduction while in contact with 
the wall. Since it is in the beginning moments of the conduction transient that most of the 
heat transfer is contributed, the initial bulk properties of the fluid are inserted in the 
expression for the transient heat flux. The familiar expression for transient conduction 
is 



1 

a 3t 


For the boundary conditions applicable to a very thick slab 


(Al) 


r = 0 T = T W 

y = 0 T = T w 

y = cc T=T M 


the solution to equation (Al) for a statistical average time period r c can be shown to be 


q 


P 


_2_ |kp£ 
Vtf T T c 


(A2) 
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and the heat -transfer coefficient h^ is 



Alternatively, equation (A3a) can be expressed as 


(A3a) 


(A3b) 


where oo is the frequency. 

The principal unknown in equation (A3a) is the contact time t c (or w). From a 
more detailed look at the mechanics of the penetration of a bulk fluid packet to the wall, 
one can gain an insight into parameters influencing the contact time r . The stages in 
the life of a packet as it comes to the wall and departs are shown in figure 6. The events 



Step 1: fresh packet 
replaces departed packet 



Step 2: fresh packet 
contacts wall 



Expanding packet \ 



Step 3: Explosive ejection 
of expanded jet 


Figure 6. - Stages in life of penetration packet. 
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shown in figure 6 begin immediately following the departure of a packet from the wall 
region. Its departure leaves a void which is immediately occupied by a new packet from 
the outer edge of the boundary layer. This is step 1 of figure 6. 

Step 2 is envisioned to be the filling of the void in the sublayer by a packet of fluid 
from the bulk region. The turbulent mixing inherent to the flow is sufficient to cause the 
void to be filled. 

The final step in the process is the rapid expansion of the fluid packet while in con- 
tact with the wall. The packet retains its identity and expands as a unit while soaking up 
heat from the wall. Only a very short contact time would be required for the entire 
packet to respond to the wall temperature environment. If the dimensions of a packet 

approximate the sublayer thickness (2x10 in. or 5x10 cm) the relaxation time for a 

6 

temperature disturbance to propagate over that thickness would be less than 10 second. 

The explosive -like expansion of the packet as it is heated from a near-critical ther- 
mal state to a temperature condition approximating wall conditions induces the inertial 
motion of fluid surrounding the packet. These inertial effects are considered to be the 
source of the jetting action that creates a void or hole in the sublayer region whereupon 
the replenishment cycle repeats again. 

From the qualitative description of the penetration mechanism presented thus far, 
the following intuitive arguments are presented in order to develop a parametric expres- 
sion for the penetration frequency. At the outset it is important to realize that the prime 
mover of the penetration model is rapid expansion by heating. Expansion amplifies a 
disturbance in the sublayer to enable a local mass exchange. It would be difficult for the 
thermal disturbance to make its influence felt in a high-speed flow. This would be ana- 
logous to the observance of the influence of a film cooling jet on a boundary layer. The 
faster the free stream traveled, the less "effective” or influenced would be a transverse 
film jet of a fixed strength. Consequently, the frequency of penetration would be inverse- 
ly proportional to the momentum in the free stream - for example, a product of the film 
density and free stream velocity. This intuition may appear to contradict the familiar 
notion that turbulence in a flowing fluid is strengthened as the flow rate (Reynolds num- 
ber) increases. In this model being proposed, I am assuming that the penetration trans- 
port is effected by an explosive -like thermal expansion of a fluid packet located in the 
sublayer. The local disruption of the sublayer can be viewed as an instability. I am 
saying that the higher shear associated with higher flow rates tends to dampen a thermal 
instability, thereby inhibiting the penetration mechanism. On the other hand, the mag- 
nitude of expansion possible within the boundary layer would be an index of the penetra- 
tion frequency. This could be portrayed by a density difference between bulk and wall 
conditions . 

In order to compare the expansion property of the boundary layer with the axial 
momentum, the two effects should be made dimensionless in a parametric form. Since 
the effects are competitive, it is convenient to cast them into a ratio form with the 
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expansion term as the numerator. Included in the expansion parameter is a theoretical 
expansion velocity, which could result from perfect conversion of thermal energy into 
kinetic energy (^j2J - H b )) • 

Since this expansion process occurs within a boundary layer, it seems appropriate 
to include viscosity terms in the parameter. The evaluations of the viscosity terms are 
made consistent with the density evaluation. Consequently, the ratio takes on the nature 
of a Reynolds number ratio: 


n . k - ?w> y^Hw - H b> (A4) 

Pf U oo Mb 

As was mentioned in the section PROPOSED HEAT-TRANSFER MODEL the previous 
dimensionless ratio modifies an adiabatic frequency w o (see eq. (4)), and can be inter- 
preted to be a frequency enhancement factor. The magnitude of cu 0 was determined 
empirically, but its value was observed to approximate measured sublayer oscillation 
frequencies. The value of u> Q used was approximately 20 cps. 

The form of equation (4) suggests that the penetration contact time is relatively long 
(low frequency) when the magnitude of the free stream momentum overpowers the expan- 
sion mechanism. A long contact time signifies an ineffective penetration mechanism or 
preservation of the sublayer. The ejection mechanism is primarily dependent on a large 
density change in the packets which penetrate to the wall. In addition, the effectiveness 
of this ejection mechanism is related to the heat flux level. This is portrayed by the 
enthalpy- difference term under the radical. Very large heating rates or extreme den- 
sity changes, or a combination of both, lead to short contact times. The shorter the con- 
tact time, the greater the frequency of the penetration mechanism, which results in an 
improved heat -transfer coefficient. 

It is interesting to observe that the amplification ratio in equation (4) is similar to 
one of the principal terms found in a stability analysis of a fluidized bed (ref. 32). The 
similarity of the penetration mechanism within the boundary layer to the mechanism of 
heat transport in a fluidized bed has been cited before in this report. Stability criteria 
for the fluidized bed as found in reference 31 are a further example wherein some simil- 
itude exists. Ruckenstein (ref. 32) suggested the pertinence of the study in reference 31 
to the penetration model of the turbulent boundary layer. Reference 32 contains a well- 
developed discussion of the penetration model and compares it to the more conventional 
time-averaged boundary layer model. 

Substituting the expression for frequency (eq. (4)) into the equation (eq. (A3a)) for 
the heat-transfer coefficient yields the useable form for heuristic comparisons with 
experimental data: 


22 




As was mentioned earlier, the reference frequency u) Q was treated as an empirical 
constant . 


Distribution of Penetration Areas 

The question does arise concerning the distribution of the areas representing the 
two heat-transfer mechanisms. In order to get some estimate of the distribution of the 
penetration areas on the surface, a simple model of the penetration mechanisms was con- 
sidered. The model assumed that the penetration mechanism was analogous to the pene- 
tration of a wall jet into a channel flow. The behavior of wall jet trajectories introduced 
normal to channel flow was studied in reference 33 and an empirical expression for the 
penetration distance was correlated. The principal dependency of the penetration dis- 
tance was with the ratio of the wall jet velocity to the free stream velocity. Using the 
Gordier correlation made it possible to develop an expression for the penetration of a 
thermal jet into the fluid stream. Perhaps the key observation in formulating the analogy 
to Gordier' s jet model is the fact that the Stanton number q/p AHU can be interpreted 
as a ratio of two orthogonal velocities. The numerator q/p AH represents a velocity 
normal to the wall and the denominator is the free stream velocity U. Such an interpre- 
tation of the Stanton number as it applies to two-phase and near-critical forced convec- 
tion is found in reference 25. For two-phase systems, it is often referred to as the 
Sterman parameter. 

The expression for penetration distance derived in appendix B and normalized with 
respect to a characteristic channel dimension is 


y u 

-E = C 9 Re 
D z 


- 0.2 

d 



(A6) 


In the interpretation of this dimensionless number, it is assumed that it represents the 
distance traveled by fluid packets as they migrate toward and away from the wall. In the 
derivation, only the motion away from the wall was considered. 

The penetration area Ap (see fig. 1) would be directly related to the penetration 
distance yp. It seems reasonable to assume that a minimal yp/D exists, which is a 
threshold of the thermally induced penetration mechanism. The penetration area fraction 
of the total surface would increase with an increase in the dimensionless penetration dis- 
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tance. A relation between the area fraction of penetration A p /A t and y /D was found 
empirically. Initially, the relation between these parameters was assumed to be of the 
form 






(A7) 


As a first guess, the exponent n was assumed to be one. It turned out through compar- 
ison of a limited amount of experimental data (11 runs) with the equation representing the 
overall model that this first approximation was a good one. Thus, the penetration area 
fraction can be represented by 



At 



Re 


- 0.2 


(A8a) 


The weak dependency on Re was dropped from the equation. Equation (A8a) was used in 
the form 


At 



(A8b) 


for values of St(L/D) < 0. 03. 

The magnitude of the coefficient was determined to be 2. 8 by using a sample set 
of hydrogen runs. Figure 7 illustrates how equation (A8b) compares to the empirical 
check of the sample data. In all of the calculations presented herein, the value of C ^ 
was fixed at the value of 2. 8. 

As is observed in figure 7, the area fraction A / A t appears to maximize at a value 
of 0. 16 and then to level off for the upper ranges of the St(L/D) parameter. This is 
somewhat analogous to what happens on boiling surfaces near burnout where the liquid- 
wetting area reaches some maximum fraction of the total area. Up to that point the 
wetting area is a function of the average heat flux (see ref. 34). 

As a further comparison of the effect of a flow rate normal to a wall interacting with 
the principal flow along a wall, some transpiration heat-transfer data were utilized. Wall 
suction by means of a transpiration surface can be considered to be a physical representa- 
tion of the penetration model presented herein. Certainly by suction, species of the free 
stream flow are carried through the boundary layer to the wall where their free stream 
molecular properties influence thermal exchange at the wall. In this respect, the suction 
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Test section, 
in. (cm) 



Figure 7. - Penetration area fraction as function of Stanton number 
and length to diameter ratio. 


process enhances heat exchange at the wall. Indeed, measurements of turbulent heat 
transfer with suction show dramatic enhancements over comparable conditions with no 
suction. This has been clearly demonstrated in the data of reference 35 and other simi- 
lar experiments. Figure 8, which is taken from reference 35, shows that the ratio of 
the transpired mass flow to the total mass flow was almost directly proportional to the 
change in Stanton number at a fixed L/D. Since the ratio of transpired flow to total flow 
can be interpreted to represent the penetration area fraction, the use of the transpiration 
data lends support to the functional relation of equation (A8b). 


Transpiration flow 
Through flow 

O No transpiration 
A 0.0024 
O .0076 



Figure 8. - Effect of transpiration flow (suction) on 
local heat transfer (from ref. 35). 
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APPENDIX B 


PENETRATION DISTANCE 

If it is assumed that the penetration distance associated with a fluid packet leaving 
the wall is analogous to the penetration distance of a jet into a fluid stream, then a rather 
simple relation for penetration distance can be developed. One experimental study of jet 
penetration normal to a flowing stream is contained in reference 33. The maximum jet 
penetration distance was observed to be a function of a ratio of the jet velocity to the free 
stream velocity. In terms of the parameters of the accompaning figure, the equation 
for maximum penetration distance is 



(Bl) 


The author of reference 33 suggested that n was 0.74. 

It is assumed that the thermal penetration distance of a packet of fluid leaving the 
wall can be correlated by the same type of equation as equation (Bl). Instead of a jet 
nozzle diameter d, the normalizing parameter is the thermal layer thickness. It is as- 
sumed that the cross-section dimension of the fluid jet emanating from the wall is of the 
order of the thermal layer thickness. Rewriting equation (Bl) as it applies to the thermal 
penetration distance and with a thermal layer thickness parameter: 

(B2) 

For simplicity n is assumed to be 0. 8. The term v/U can be evaluated in several 
ways. One way is to recognize that v/U = f(St) = q/AHpU where AH is enthalpy differ- 
ence between wall and bulk condition. Modifying equation (B2) to include channel height 
yields 


^max 

D 



(p (gt) 0. 8 


(B3) 
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For a constant q surface, a simplified treatment of the integral energy equation in 
axisymmetric coordinates yields the energy thickness (information obtained from a pri- 
vate communication with D. R. Boldman of Lewis) as follows: 

£ = 0.031 Re^ 0 - 2 Q^°’ 8 (B4) 

Equation (B4) depicts the development of the thermal layer in the entrance region along 
a flat plate or along a tube. Experimental studies have shown that the thermal entrance 
effects may persist as far down as L/D = 40. Substitution of equation (B3) into (B4) 
yields 

S (St)°- 8 (B5) 

This equation is assumed to retain the L/D correction up to values of L/D - 40. There- 
after, the L/D correction is held constant at the L/D - 40 value. 
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